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ABSTRACT 

In modeling the rapid part of the pressure correlation term in the Reynolds stress 
transport equations, extensive use has been made of its exact properties which were first 
suggested by Rotta 1 . These, for example, have been employed in obtaining the widely 
used LRR model. Some recent proposals have dropped one of these properties to obtain 
new models. We demonstrate, by computing some simple homogeneous flows, that doing 
so does not lead to any significant improvements over the LRR model and it is not the 
right direction in improving the performance of existing models. The reason for this, in our 
opinion, is that violation of one of the exact properties can not bring in any new physics 
into the model. We compute thirteen homogeneous flows using LRR (with a recalibrated 
rapid term constant), IP and SSG models. The flows computed include the flow through 
axisymmetric contraction; axisymmetric expansion; distortion by plane strain; and homo- 
geneous shear flows with and without rotation. Results show that for most cases a properly 
calibrated LRR model, which is the most general representation for a model linear m the 
anisotropic tensor, performs either better or as good as the other two models of the same 


level. 



1. INTRODUCTION 


Modeling of the pressure strain correlation is one of the central issues in the develop- 
ment of Reynolds stress closure models. Most of the popularly used models have utilized 
the suggestions of Rotta 2 in developing models for this correlation. From the Poisson 
equation for the fluctuating pressure he obtained the following expression for the pressure 
correlation 
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where the superscript S refers to the so called slow part and the superscript R refers to 
the rapid part of the pressure correlation 2,3 ’ 4 . If we restrict ourselves to the homogeneous 
turbulence with constant mean gradients (only such flows are considered in this study) 
then U p , q can be moved outside the integral and one has to model only what else is left 
inside the integral. For brevity we designate 
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Rotta showed that the fourth rank tensor I qip j has the following properties which are exact. 


Iqipj = I iqpj = Iqijp, symmetry property 
Iqipi = 0, incompressibility, 

Iqijj = u 9 Uj, normalization property (4) 

Launder et al 2 utilized these constraints in obtaining the so called Launder, Reece and 
Rodi (LRR) model. They first assumed the most general form of Iq ipj which is linear in the 
anisotropic tensor b {i = (uiuj/q 2 - S^/Z). Then application of the symmetry constraint 
reduced the number of unknown constants to five. Further application of the incompress- 
ibility and normalization constraints reduced the number of unknown constants to one 
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which must be determined empirically. Recent study of Speziale et al. 5 has proposed a 
model of Iqipj which is linear in the anisotropic tensor bij and satisfies all of the above con- 
straints except the normalization constraint. This model, therefore, involves four empirical 
constants. If the normalization constraint is applied to the SSG model then it reduces to 
the LRR model. So the most general rapid term model which is linear in the anisotropic 
tensor and satisfies the above mentioned constraints is the LRR model. In our opinion 
the violation of the normalization property of the rapid term correlation does not bring 
in any new physics to the model and, therefore, is not the right direction in improving 
the performance of the currently used models. Direct comparisons of rapid term model 
with DNS data 6 showed that the LRR model performs better than SSG model. In this 
paper it will be further demonstrated, through computations of homogeneous flows, that 
the properly calibrated LRR model can perform as good as SSG and IP models or even 
better. 

The thirteen homogeneous flows that are numerically computed include flow through 
axisymmetric contractions and expansions; distortion by plane strain; homogenous shear 
flows with and without rotation of coordinate axis. In the light of the objective of this 
study we concentrate only on those models which use a linear (or quasi-linear) model for 
the rapid term. This, therefore, includes the Launder , Reece and Rodi 2 (LRR) model, 
the isotropization of production (IP) model 2 , and the Speziale, Sarkar and Gatski 5 (SSG) 
model. In the next section governing equations and the models tested are given for ho- 
mogeneous flows. Then a brief description is given of the flows computed followed by the 
results and conclusions of the study. 

2. GOVERNING EQUATIONS AND MODELS 

For an incompressible homogeneous flow in a rotating frame the Reynolds stress equa- 
tion can be written as 

— tlTEJ = - (ujUfc Uj t k + ujul U it k ) - 2(v^u^e m kjCl m + uJuTe m fe»O m ) 

+ ^P(«»,j + «i,i) - 2 vvT^UjJ (3) 
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where fl m is the rotation rate vector of the coordinates relative to the inertial frame. The 
model for the pressure strain correlation (denoted by II »j) and a modeled dissipation rate 
tensor are needed to complete the second order closure. Usually, the deviatoric part of 
dissipation tensor is combined with the pressure strain, and a dissipation rate (i.e. e) 
equation is needed. 

The form of the dissipation rate equation used by all the models tested here is given 
by 

Dc e e 2 

- Cti-UiUkUi'k - C t -z— (4) 

where k = ujui/ 2; C ( \ and C ( 2 are model coefficients. For LRR and IP models we use 
C ( 1 = 1.44 and C ( 2 = 1-92 whereas for SSG model C ( 1 = 1.44 and C f2 = 1.83. 

LRR MODEL 

The LRR model for the pressure correlation term is given by 

Hij = — Citbij 

+ 0.8 kSij + ~k(b ik Sjk + bjkSik - | b mn Smnf>ij ) 

+ (2 ° + %>**) (5) 

where b {j = (uiuj/q 2 - Sij/ 3), 5 i; - = (U itj + U jti )/2 and W {j = (U itj - U jfi )/2 + e mji Sl m . 
The model constants of the LRR model have evolved to slightly different values than those 
originally recommended by LRR 2 . The value of the Rotta constant Ci (in the return to 
isotropy term) used in the present study is 3.6 This value was used, among others, by 
Gibson and Launder"^. The value of the rapid term constant C\ (is 0.4 in the original 
LRR model) used in this study is 0.55, which is slightly higher than the value of 0.5 
recommended by Morris 8 . Therefore, in this paper LRR model will mean equation (5) 
with C\ = 3.6 and C 2 — 0.55. It was found that the value of 0.55 led to improvement in 
the performance of LRR model for all the flows tested here. (The improvements were slight 
for the irrotationally strained flows but were significant for the homogeneous shear flows.) 
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It should be mentioned that originally LRR obtained the value of this constant from the 
experiments of Champagne et al 9 which were later redone under higher shear rate (see e.g. 
Tavoularis and Corrsin 10 ) to obtain more consistent results and these later experiments 
have been used by several researchers 11 as a standard benchmark for calibrating model 
constants. A slightly different value of rapid constant used here is merely a reflection of 
this fact. 

The LRR model is linear in Reynolds stress and satisfies the constraints of incom- 
pressibility, normalization and symmetry. It does not satisfy the realizability principle ( 
which implies that the non-negative quantities such as u 2 should remain non-negative, see 
Shih and Lumley 4 for details). It should be further noted that the LRR model is the most 
general model for the rapid term which is tensorially linear in the anisotropic tensor and 
satisfies the above mentioned three basic constraints. 

IP MODEL 

This model uses the same slow term model as the LRR but for the rapid term keeps 
only the — 2P6ij/3) term of the LRR model. When re-written in terms of Wy, 

and Sij this model is given by the following expression. 

2 

lit; = — Cjtbij + C 2 ^k [(fcjJtS'jJt + bjkSik -bmnSmnfiij) 

+ (b ik W jk + b jk W ik ) + p ij ] (6) 

where C\ = 3.6 and C 2 = 0.6 were recommended by Gibson and Launder 7 . 

SSG MODEL 

The SSG model for the pressure correlation term is given by 

n a = - (<7,e + c;p)bij + c 2 e(b ik b kj - + (c 3 - c;ii'' 2 )kSij 

+ C4k(bi k Sjk + bjkSik - -b mnSmubij C 5 k(b ik w jk + b jk W ik ) (7) 

where the definitions of bij , S{j and W x j are the same as given earlier; P = —uiujUij and 
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II — bijbij. This model involves seven empirical constants and their values are given by 

Cj = 3.4, Cj* = 1.80, Ci — 4.2 
C 3 = |,<7 3 * = 1.30,(74 - 1.25, C 5 = 0.40 

The rapid part of the SSG model satisfies the constraints of incompressibility and symmetry 
but it does not satisfy the normalization property. If the constraint of normalization was 
imposed, the rapid part of SSG model would reduce to LRR model. 


FLOWS COMPUTED 

A total of thirteen flows are computed in the present study. These include two ax- 
isymmetric contractions (DNS) 12 ; two axisymmetrfc expansions (DNS) 12 ; four distortions 
by plane strain (DNS) 12 ; two homogeneous shear flows (one DNS 13 and one experiment 10 ); 
and three homogeneous shear flows with rotation (LES) 14 . For the flow through the ax- 
isymmetric contraction the mean velocity gradient is given by 



0 0 \ 

~im o 1 

2d* U I 

0 -I»s / 

2 Bx / 


w ^ere T* 1s positive. For the flow through axisymmetric expansion the mean velocity 
gradient tensor is the same as above except that ^ is negative. For distortion by plane 
strain flow the mean velocity gradient tensor is given by 



where is negative. For the homogeneous shear flow the mean velocity gradient tensor 
is given by 


Uij = 


0 f 0 
0 0 0 


0 0 0 

For the homogeneous shear flow in a rotating frame the mean velocity gradient tensor is 
same as above but in addition it has a non-zero rotation rate which for the cases considered 
here is given by Dj = (0, 0, D) where Q is constant and represents the rotation rate of the 
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non-inertial frame of reference relative to the inertial frame of reference. In all these cases 
the mean velocity gradients are constant. All these flows constitute an initial value problem 
and were numerically solved using a fourth order Runge-Kutta finite difference method. 

3. RESULTS 

Figure 1 compares the development of Reynolds stresses computed using these three 
models in a flow through axisymmetric contraction with the DNS data. (Following Lee 
et al. 12 all the stresses will be plotted as a function of the normalized time St.) Here we 
show a typical case of S = 10.00 ( Sk 0 /e 0 = 5. 573, case AXL). All the models deviate from 
the DNS data. However, LRR model gives slightly better results than the SSG model for 
both the components while IP model performs the worst. Figure 2 shows a comparison 
for a much higher strain rate case ( S = 100.00, Sk 0 /e 0 = 55.73, case AXM). We note that 
now all the models deviate even more from the DNS data. We also observe that the LRR 
model gives better results than the SSG model. 

Figures 3 and 4 show a similar comparison for flow through axisymmetric expansion for 
two different strain rates 12 . For the smaller strain rate flow ( S = 0.717, Sk 0 / e Q = .408, case 
EXO) SSG model reproduces the u 2 development quite well while both IP and LRR models 
under-predict it. For the v 2 component all the models give similar results. Therefore, for 
this low strain rate flow SSG model is better than the other two models. For the flow 
with higher strain rate ( S = 7.17, Sk 0 /e 0 = 4.08, case EXP) the LRR model is in excellent 
agreement with the DNS data for both the components while both IP and SSG models 
show over-prediction So for this flow LRR model works the best. 

Now we show comparisons for the distortion of turbulence by plane strain for four 
cases of differing strain rates 12 . We start from the smaller strain rate case. Figure 5 
compares the evolution of the three non-zero Reynolds stress components for the flow with 
strain rate 5 = 0.65 ( Sk 0 /e 0 = 0.577, case PXA). We observe that all the models predict 
u 2 component reasonably with LRR model slightly better than the other two models. For 
the v 2 component SSG model is closest to the data where as both LRR and IP models 
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slightly under-predict it. For the w 2 component all the three models give about the similar 
results and predict it reasonably. Figure 6 shows the comparison for the case of S — 1.3 
( Sk 0 J e 0 = 1.154, case PXB). The u 2 and w 2 components are reproduced better by the LRR 
model. However, for the v 2 component SSG model gives better prediction than LRR and 
IP models. The comparisons for the third case, with strain rate S = 2.6 ( Sk 0 /e 0 = 2.309, 
case PXC), are shown in figure 7. For u 2 component all the models under-predict the DNS 
data. LRR model is slightly better than the SSG model. IP model is the worst of the 
three. For v 2 component IP model works the best. LRR model slightly under- predicts v 2 
while SSG over-predicts it. The third component w 2 is over-predicted by all the models 
with LRR model being better than the other two. Figure 8 shows the similar comparisons 
for the highest strain rate case ( S = 25.0 ,S 0 k/e 0 = 22.227, case PXE). All the three 
models under-predict the u 2 component. IP model is the worst of the three models. LRR 
model gives slightly better result than the SSG model for this stress component. For v 2 
component LRR model is the best and SSG model is the worst of the three. For the w 2 
component all the three models over-predict the DNS data with LRR model being closest 
to the data. From the above four plane strain flow comparisons, we note that on the overall 
LRR model works better than the other two models. 

Now we turn our attention to homogeneous shear flow. First we show the comparisons 
with the DNS data 13 for the shear rate S = 56.56 ( Sk 0 /e 0 = 2.36, case C128W). Figure 9 
shows that both LRR and SSG models preform equally quite well. On the other hand IP 
model is off the data for all the components. 

Figure 10 shows the same comparison with the homogeneous shear flow experiment 
of Tavoularis and Corrsin 10 (5 1 = 46.8, Sk 0 /e 0 — 6.46). (Following Tavoularis and Corrsin 
all the Reynolds stresses shown here are normalized by mean velocity along the centerline 
of their wind tunnel and the streamwise distance is normalized by the height of the wind 
tunnel.) For the u 2 component LRR model gives the best result whereas SSG and IP 
models over-predict it. For the v 2 component also the LRR works the best. SSG model 
slightly over-predicts the data whereas IP model is off by a larger margin. For the w 2 
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component both SSG and IP models reproduce the data very well whereas LRR model 
over-predicts the data. For the shear stress component LRR performs very well whereas 
SSG model over-predicts the data and IP model is off the data by a larger margin. So 
for this experiment LRR model has better overall performance than the other two models. 
Finally, we discuss the evolution of q 2 for the case of rotating homogeneous shear flow. 
Since no experimental or DNS data is available for this flow the comparisons will be made 
(for two cases) with the LES of Bardina et al. 14 . Bardina 15 has pointed out that in this 
case we should be careful in interpreting the comparisons and should use the LES data 
for judging the model performance in a qualitative fashion. In all the cases shown here 
the initial conditions corresponded to isotropic turbulence with e 0 /Sk o = 0.296. Figure 
11 shows the comparisons for the three cases of different Rosby numbers (= Cl/S). For 
Cl/S = .25 we note that all the three models significantly under-predict the LES results 
for q 2 ; SSG being closest to the LES data and the LRR being the furthest. Qualitatively 
all the three models reproduce the LES trends. For the case of Cl/S — 0.50 SSG is in 
excellent agreement with the LES results. Both IP and LRR give identical results and give 
a smaller value of q 2 than the LES. For the third case of Cl/S = 1.0, all the three models 
give identical results. Since no LES results are available for this case the only purpose of 
showing the results is to see how the three models compare with each other. 

4. CONCLUSIONS 

Results were shown from numerical computation of various homogeneous turbulent 
flows using three different turbulence models. All of these models use a linear (or quasi- 
linear) model for the rapid part of the pressure strain term. The rapid term constant in the 
LRR model was recalibrated to a slightly different value then the one originally proposed 
by LRR. Based on their overall performance it is found that the recalibrated LRR model 
works as good as the SSG and IP models. For the irrotational flows the differences the 
performance of the LRR model improved over the other two models as the strain rate 
increased. For the simple homogeneous shear flow experiment LRR model worked better 
than the SSG model. For the DNS of the same flow both performed equally good. For 
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the rotating homogeneous shear flows both SSG and LRR model showed trends similar 
to those shown by LES with SSG performing closer to LES. It is worth noting that SSG 
model has seven empirical constants as compared to two in LRR model and still, on the 
overall, it does not perform better than LRR model. Part of the reason for this may be 
due to the fact that the SSG model does not satisfy the normalization constraint where 
as LRR model does. For a model of the rapid pressure strain term which is tensorially 
linear in the anisotropic tensor and satisfies all the three basic constraints LRR is the most 
general model. The comparisons shown in this study demonstrate that the model obtained 
by relaxing normalization property of the pressure correlation term does not lead to any 
significant improvement over the LRR model and, therefore, this is not the right direction 
in improving the performance of the current models. 
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Figure 1. Comparison of the models for flow through axisymmetric contraction with the 
DNS data of of Lee and Reynolds (1985), case AXL, S=10.0. 
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Figure 6 (continued from last page) 
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